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The stellar (n, y) cross section data for the mass numbers around A œ% 160 are of key importance 
to nucleosynthesis in the main component of the slow neutron capture process, which occurs in the 
thermally pulsing asymptotic giant branch (TP-AGB). The new measurement of (n, y) cross sections 
for P9'Tb was performed using the CgDg detector system at the back streaming white neutron 
beam line (Back-n) of the China spallation neutron source (CSNS) with neutron energies ranging 
from 1 eV to 1 MeV. Experimental resonance capture kernels are reported up to 1.2 keV neutron 
energy with this capture measurement. Maxwellian-averaged cross sections (MACS) are derived 
from the measured P?'Tb (n, y) cross sections at kT = 5-100 keV and are in good agreement with 
the recommended data of KADoNiS-v0.3 and JEFF-3.3, while KADoNiS-v1.0 and ENDF-VIII.O 
significantly overestimate the present MACS up to 40% and 20%, respectively. A sensitive test of 
the s-process nucleosynthesis is also performed with the stellar evolution code MESA. Significant 
changes in abundances around A ~ 160 are observed between the ENDF/B-VIII.0 and present 
measured rate of 1°°Tb(n, y)'®°Tb in the MESA simulation. 
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I. INTRODUCTION 


The elements heavier than iron in the solar system are mainly produced by the slow neutron capture process (s- 
process) [1] and the rapid neutron capture process (r-process) [2] in stars as found by Burbidge et al. [3] and Cameron [4] 
as early as in 1957. Almost less than 196 of heavy elements are ascribed to the production of photo-disintegration 
processes (the so-called p/y process), neutrinos, and charged-particle induced reactions [5, 6]. The s-process takes 
place during stellar evolution and path through nuclei along the valley of 8 stability with lower neutron densities and 
temperatures. The weak s-process is responsible for producing isotopes up to A ~ 90 and occurs during He and C 
burning in massive stars. While the main s-process contributes isotopes A ~ 90-208 and take place in the thermally 
pulsing asymptotic giant branch (TP-AGB) phase of low- and intermediate-mass stars. The r-process is related to 
explosive nucleosynthesis in massive stars and binary star mergers with relatively high neutron densities. 

Terbium is mainly produced by the explosive r-process [2], while about 9% is made by main s-process in low- and 
intermediate-mass TP-AGB stars [7]. The s-process reaction path around terbium is sketched in Fig. 1. 1©°Tb is an 
important s-process branching point, which is shielded against the G~ decay chains from the r-process by stable nuclei 
160Gd. The abundance of !9"Tb will affect the s-only isotope !99? Dy, which paths through the s-process to !6! Dy and 
considered as an reference point for the reaction flow. Also, neutron capture of terbium isotopes will contribute to the 
production of p-nuclei ‘Er through the temperature sensitive branching at !6?Dy to 163Ho by the nuclear reaction 
path 1©3Ho(n,y)'64Ho(8)!®Er in the stellar environment [8-10]. 

Two groups of experimental neutron capture cross section data for !99'Tb(n, y) reaction were reported between 1961 
and 1991 [11-16] in the astrophysical interested energy region. Data from Block et al. [11], Gibbons et al. [12], Lépine 
et al. [13] are significantly larger than data of Mizumoto et al. [14], Ohkubo et al. [15] and Bokhovko et al. [16]. In the 
KADoNiS-v0.3 [17] database, the recommended Maxwellian-averaged cross sections (MACS) of °°Tb(n, y) at stellar 
temperature kT = 30 keV is 1580 + 150 mb, which is widely used in stellar nucleosynthesis calculations, derived from 
the neutron capture cross section data of Mizumoto et al. [14] and Bokhovko et al. [16]. While, this value reported by 
the recent activation measurement [18] with 2166 + 181 mb and it varies from 1116 to 2949 mb in the collection list 
of the KADoNiS-v1.0 [19] and references therein. Because of so large discrepancies, newer measurements of stellar 
neutron capture cross section of ?'Tb are desired for improving the reliability of recommended data. Much higher 
accurate neutron capture cross sections for the terbium isotopes are required to the s-process nucleosynthesis, and the 
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FIG. 1: (color online) The s-process reaction path in the region of terbium. Black solid boxes indicate stable; blue empty boxes indicate 
unstable isotopes. Arrows to the right and the next higher elements represent neutron capture reactions and f^ decays, respectively. The 
isotope! Dy is a terrestrially stable nucleus, which becomes unstable at stellar temperature [8]. 


separation of the r and s components in the observed Tb abundance through subtraction of the s abundances from 
the respective solar values. In addition, the !99?'Tb(n, y) cross section data is very important for the design of nuclear 
reactor and the research in nuclear structure [20, 21]. 

In the past years, a series of neutron radiative capture experiments have been performed using the CgDg detectors 
with the total energy technique [22-26] and the total absorption calorimeter consists of BaF» scintllator arrays [27— 
30] in combination with the "Li(p, n)'Be reaction neutrons and white neutron sources. Both of the CeDe and BaF» 
detection systems, respectively, were constructed at the back streaming white neutron facility (Back-n) of the China 
Spallation Neutron Source(CSNS) [31, 32], which is the only high intense pulsed spallation reaction neutron source 
in China at present. The main physics motivation of these systems is to measure the neutron capture data related 
to the nuclear astrophysics, advanced nuclear energy technologies such as accelerator-driven subcritical systems, 
thorium-based molten salt reactor, and fourth generation reactors. The detection systems, experimental techniques, 
background study and resonance energy region measurements have been described in Refs [33-36]. In this work, we 
report the new experiment of neutron capture cross sections for 159Tb over the energy range from 1 eV to 1 MeV via 
the time-of-flight (TOF) method combined with the white neutron beam at the CSNS Back-n facility. The resonance 
kernels of °°Tb (n, y) are obtained up to 1.2 keV, and the MACS for the astrophysical interest energy region from 
kT = 5-100 keV are calculated. 


II. EXPERIMENT 


The measurement of the !9?'Tb (n, y) was performed at the CSNS Back-n facility [31, 32] in China. Neutrons were 
produced via the spallation reactions by bombarding a massive tungsten target with a 1.6 GeV proton beam, operating 
at a 25 Hz repetition rate and collimated to a 30 mm diameter beam spot at the sample position. The neutron energies 
ranging from thermal to a few hundred MeV were determined using the time of flight ((TOF) method with a 76 m 
long flight path. The neutron capture yields were studied using the total energy measurement principle based on an 
array of four CeDe detectors with the pulse height weighting techniques (PHWT) [37-39]. The neutron flux in the 
energy range of 1 eV to 100 MeV was obtained with a combination of the data from dedicated measurements of the 
°LiF-silicon (Li-Si) detector array [40] and the ?35U loaded multilayer fission chamber [41]. The Li-Si detector neutron 
flux was normalized to ???U data at 10-20 keV since its position was about 20 m away from the sample. The samples 
used for this experiment are listed in Table I. A metallic terbium sample of natural abundance with a thickness of 
0.2 mm was used for the 15°Tb (n, y) reaction study. The !?" Au sample was used as the standard for calculating the 
relative cross section of 159Tb (n, y) in neutron energies ranging from 2 keV to 1 MeV. The "Pb samples were used 
to normalize the neutron fluence and to evaluate the in-beam y-rays and scattered-neutron backgrounds, respectively. 
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'The empty sample runs were also performed to evaluate the background produced from the upstream devices, such 
as the sample holders, etc. All experimental data were collected on an event-by-event basis high-performance data 
acquisition system with a 12-bit full-waveform digitizers sampling at 1 GS/s [42]. The offline data analysis was done 
on the CERN ROOT framework [43]. 


TABLE I: Characteristic parameters of samples 


Sample Thickness (mm) Diameter (mm) Mass (mg) Area desity (atom - b^!) 
nei 0.20 30 1169.14 6.27 x 1074 
197 Au 0.10 30 1357.17 5.87 x 107* 
EP 0.53 30 4249.75 1.75 x 107° 


Empty holder 


III. DATA ANALYSIS AND RESULTS 


The incident neutron energy £j was determined by employing the TOF method through the nonrelativistic kine- 
matics formula (1). The neutron effective flight path L was obtained by analyzing the low energy resonances of the 
gold sample. 


72.2977 x L 
En = (MÁY (1) 
where E, is in MeV, L is in meters, and t, is the flight time in ns. 
The experimental neutron capture yield as a function of E, can be calculated as: 


_ 1C"(E.) - BY (En) 


Yezp(En) = f D(En) x E, (2) 


where f is the normalization factor determined by self-normalizing the measured capture yield of 4.9 eV resonance 
of 1?" Au and 11.1 eV resonance of !159Tb, based on the saturated resonance technique [44, 45]. E. is the detection 
efficiency for a capture event, ®(£,,) is the neutron flux spectrum, and C'"(E,) and B” (En) are the weighted count 
spectrum and the total background spectrum, respectively. 

The PHWT method is essential for the (n,y) cross section measurement using the CgDg detection system. Its 
function is to make the detection efficiency £4 proportional to the incident y ray energy Æ}, as Eq. (3) 


Ey = L W F;(Ea)R(Ea, Ey) = aE, (3) 


where a is the constant parameter, R is the detector response function, E; is the energy bin of pulse height spectrum, 
and WF is the weighting function, which can be approximated by a 5'^" polynomial function. The polynomial 
coefficients are obtained with a minimum least-squares fit to the the detector response for 27 different monoenergetic 
y rays from 0.1 MeV to 10 MeV simulated by the GEANT4 code [46], taking into account the detailed model of the 
detection systems and the sample. 

The background consists of two components, sample-independent and sample-dependent ones. The former is eval- 
uated by the empty sample runs under the same experimental conditions. The sample-dependent background is 
mainly induced by in-beam ^ rays and scattered neutrons on the sample. Their contributions were determined by 
measurements with a "**Pb sample and parameterized as eq.(4). 


where the B, and B, denote the background contribution of the in-beam gamma rays and scattered neutrons, 
respectively, and can be formulated as eqs. (5-6). 


B (En) =b x e-*/ VE" +d x ec eX VE 4 f (5) 


(6) 


In order to normalize these background components, By and Bn, dedicated runs with the !5! Ta and ??Co neutron 
filters were performed. Fig. 2 shows the energy spectrum obtained from the °Tb with filters and background 
components. The normalization factors f, and fy for Bn and B, components were obtained by matching the dips 
of the filtered spectra. The effect of the filters on in-beam gamma rays and neutrons is evaluated by considering the 
neutron flux and energy distribution of in-beam gamma rays. The neutron and gamma energy spectra of Back-n are 
sampled for the incident particle energy spectra using GEANTA code to simulate with and without filters, and the 
counts of scattered neutrons and gamma are recorded at the detector position [33]. The reduced neutron and gamma 
attenuation factors are 0.92 and 0.68, respectively, and are used for corrections to fn and fy. 
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FIG. 2: (color online) 9'Tb(n, y) capture yield with filters and background components are shown as a function of neutron energy. 
Black and blue solid lines indicate the !99?'Tb sample with 19! Ta + 9? Co neutron filters and total background which is made of scattered 
neutrons (green dash-dotted lines) and in-beam y rays (red dashed lines). 


The resonance parameters (i.e., the radiative width T, and the neutron width I) are extracted using the R-matrix 
code SAMMY [47] to fit the experimental capture yields of °°Tb(n, y) in the resonance region up to 1.2 keV with 
the JEFF-3.3 [48] evaluated data as initial values. The code includes the experimental conditions, such as sample 
composition, multiple scattering, neutron self-shielding, and the Doppler effect at room temperature. The energy- 
dependent parameterized resolution function of the Back-n facility [49] is used to broaden the resonance peaks. The 
measured capture yields together with the SAMMY fits in neutron energy ranges below 100 eV are shown in Fig. 3. 
Good agreements were observed between present data and SAMMY fits both in the resonance energy and shape. The 
resonance parameters obtained in this work are compared with those from JEFF-3.3 evaluations given in Appendix 
A. The experimental capture resonance parameters show substantial agreement with the JEFF-3.3 evaluations in the 
energy range below 100 eV. However, differences are noted in the energy range between 100 eV and 1.2 keV, which 
is attributed to the degradation of the experimental resolution of the Back-n facility in this measurement. Average 
resonance parameters (such as the average level spacing, the neutron strength functions, the average radiation widths) 
are determined by fitting the measured average capture cross section data in the energy range between 2 keV and 200 
keV with the code FITACS [50], implemented in SAMMY. The resonance parameters given by Mughabghab [51] are 
used for initial values in the fitting procedure and the average level spacing is specified as 3.78 eV [51]. The results 
obtained from the FITACS fit and the recommended data [51, 52] are given in IX. 
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The measured neutron capture cross section ratio of !^9?Tb relative to the standard !?"Au sample is obtained in 
the neutron energy range of 2 keV to 1 MeV using the same method of Refs. [29, 53]. The absolute cross sections 
are converted from these ratios using the gold data of JEFF-3.3 evaluations. In Fig. 4, the experimental !9?'Tb(n, 
y) cross section and TALYS-1.9 [54] calculation result of this work are compared with the available measured data 
in literature [11-16] and evaluated data in libraries JEFF-3.3 [48], ENDF/B-VIIL.0 [52]. As illustrated, the present 
experimental results are in good agreement with the data from TALYS-1.9, Bokhovko et al. [16], and Ohkubo et 
al. [15] in the overlapping energy regions. The cross section of JEFF-3.3 and Mizumoto et al. [14] are best fitted by 
the present data below 400 keV neutron energy and slightly larger above 400 keV. On the other hand, the data of 
Block et al. [11], Gibbons et al. [12], Lépine et al. [13] and ENDF/B-VIII.0 library are significantly overestimated 
result of this work. 
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FIG. 3: (color online) The SAMMY fits to the experimental capture yields of 1°°Tb (n, ^). 


The Maxwellian-averaged cross sections (MACS) for stellar temperature can be calculated from the neutron capture 
cross sections using the formula (7) [23, 30]: 


2 sa -Ejk 
MACS = TOTE f c(E)Ee-F/*T dE (7) 


where kT is stellar temperature, o (E) is the neutron capture cross section, E is the neutron energy. 


TABLE II: The contributions from JEFF-3.3 evaluated data to present MACS. 


kT (keV) 5 10 20 30 40 50 60 TO 80 100 
Ratio (%) 10.0 4.0 1.5 0.85 0.56 0.41 0.32 0.28 0.21 0.15 


For °°Tb(n, y), the present MACS are determined using three parts of cross sections from 1 eV to 1 MeV energy 
range. The first part, shown in Fig. 5(a), is the cross sections reconstructed from the SAMMY fits to the present 
measured neutron capture yield in the resolved resonance region below 100 eV. The second part, shown in Fig. 5(b), 
is derived from the JEFF-3.3 evaluated cross section data between 100 eV and 2 keV since the resonance parameters 
obtained in this study are not sufficiently reliable due to the degradation of the experimental resolution. The JEFF-3.3 
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FIG. 4: (color online) Comparisons of the measured neutron capture cross section on !?9'Tb of this work (black full circles) with previous 
measurements [11-16] (symbols), evaluations [48, 52] (dashed lines) and TALYS-1.9 calculations [54] in the 2 keV to 1 MeV energy region. 


data contributions to the whole MACS are found in Table II. As we can see, this component cannot be ignored in 
the current study of s-process nucleosynthesis. The last part is provided by the present experimental cross sections 
relative to the gold standard over the energy range of 2 keV to 1 MeV as shown in Fig. 5(c). 

The result of this work at kT = 30 keV is compared with previous experiments, theoretical calculations, and 
recommended data from KADoNis database in Table III. The MACS recommended data from KADoNiS-v0.3 and 
the calculations of TALYS-1.9 code agree with the present measured data less than 496. The values from the recent 
activation measurement of Praena et al. [18], and the earlier data of Allen et al. [55], Lépine et al. [13] are relatively 
larger than the present result, whereas the results of Bokhovko et al. [16] and NON-SMOKER [58] are smaller. 
Theoretical calculations from MOST 2002 [57] and MOST 2005 [56] are obviously overestimated and underestimated 
by 8996 and 2896, respectively, compared to the present data. 
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FIG. 5: (color online) The °° Tb(n, y) cross sections used in the calculation of present MACS. 


The MACS values are better to be known as a function of thermal energies between a few to hundred keV at the 


202306.00617v1 


chinaXiv 


ẹ This Work 

th KADoNis-v0.3 

o KADONis-v1.0 
— TALYS-1.9 
[ — — ENDF/B-VIII.0 
---- JEFF-3.3 


10^ 


0 20 40 60 80 100 
kT (keV) 


FIG. 6: (color online) (top) The MACS of 159Tb(n, y) and (bottom) the ratios between present experimental data and others [17, 19, 48, 


52]. 


TABLE III: The present MACS (in mb) of ?Tb(n, y) at kT = 30 keV compared to past experiments, calculations and 
evaluations. 


This work 1560+160 KADoNis-v0.3 [17] 1580+150 KADOoNiS-v1.0 [19] 1817+258 
Praena 2014 [18] 2166+181 Bokhovko 1992 [16] 1471+66 Lépine 1972 [13] 1850+250 
Allen 1971 [55] 2200+200 MOST 2005 [56] 1116 MOST 2002 [57] 2949 
NON-SMOKER [58] 1427 TALYS-1.9 [54] 1628 


TABLE IV: The derived MACS for '°?Tb(n, y) in mb. 
kT(keV) This work KADoNis-v0.3 [17] KADoNiS-v1.0 [19] TALYS-1.9 [54] ENDF/B-VIILO [52] JEFF-3.3 [48] 


5 41314421 4365 48754851 3892 5726 4028 
10 2814+287 2930 3323573 2790 3893 2750 
15 2262+231 2330 2664+431 2294 3092 2232 
20 1937198 2010 2276+348 2014 2623 1927 
25 1715175 1750 2012+295 1792 2307 1716 
30 1551158 1580+150 1817+258 1635 2076 1559 
40 1318+134 1340 1544211 1395 1756 1334 
50 1159118 1185 1358181 1225 1540 1178 
60 1041106 1070 12204159 1098 1381 1063 
80 877+91 910 1029+128 920 1158 902 
100 766482 800 898+106 811 1005 794 


various s-process sites. For completeness, the comparisons of the experimental MACS and the TALYS-1.9 calculation 
result of this work to the recommended data of KADoNis and evaluation libraries for stellar thermal energy from 
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5 to 100 keV are shown in Fig. 6 and Table IV. On top panel: the black full circles and red solid line denote the 
present measured data and TALYS-1.9 calculation, respectively; the open green squares and blue circles indicate 
KADoNis-v0.3 and KADoNis-v1.0 recommended data; the blue dashed line and pink dotted line refer to calculated 
MACS by using the evaluated cross section data of ENDF/B-VIILO and JEFF-3.3 libraries. On bottom panel: the 
ratios of the measured MACS to those from KADoNis-v0.3 (green open squares), KADoNis-v1.0 (blue open circles), 
ENDF/B-VIII.0 (blue dashed line), TALYS-1.9 (red solid line), and JEFF-3.3 (pink dotted line) data; Black solid line 
corresponds to ratio 1. Our measured MACS data is well reproduced by the TALYS-1.9 calculation and the JEFF-3.3 
evaluation less than 596 for the entire thermal energies studied. While, the present MACS is significantly smaller than 
those of ENDF/B-VIIL.0 evaluation by 3096—4096. It should be noted that, the present experimental MACS is smaller 
than KADoNiS-v0.3 data by about 396-1196 and 4% for the energy region kT < 20 keV and 20 keV < kT < 100 keV, 
respectively. MACS data derived from KADoNis-v1.0 is larger about 1596-2596 than the values of this work in the 
whole energy range. 

The total experimental uncertainties of MACS include statistical and systematic contributions as shown in Table V. 
The statistical uncertainty is less than 1% below 100 eV and less than 5% between 2 keV and 1 MeV. The systematic 
uncertainties are due to neutron flux (4.596 below 150 keV and 8.096 above 150 keV), saturated resonance peak 
normalization (196), background subtraction with filters (8.696), flight path (0.0896, neglected), PHWT calculation 
(396) and sample impurity (0.0196). Additionally, the contributions from neutron inelastic to the present MACS are 
included at different stellar temperatures as shown in Table VI, which cannot be eliminated in the present experiment. 


TABLE V: Systematic uncertainties of MACS ??Tb (n, y) cross sections. 


Source Uncertainty (96) 
Neutron flux («150 keV; 2150 keV) 4.5; 8.0 

Pulse height weighting functions 3.0 
Background subtraction with filters 8.6 
Normalization factor 1.0 

flight path 0.08 

Sample impurities 0.01 

Total (without inelastic contribution) 10.2;12.2 


TABLE VI: The contributions from neutron inelastic to present MACS. 


kT (keV) 5 10 20 30 40 50 60 70 80 100 
Ratio (%) 107°} 10-! 107° 107% 107? 0.21 0.68 1.61 3.11 7.90 


The astrophysical reaction rates at a given thermal temperature can be numerically calculated from corresponding 
MACS as eq.(8): 


Na < ov > (kT) = Na x MACS xv (8) 


where N4 is the Avogadro number, MACS is defined by Equation (7), v = \/2kT/ is the mean thermal velocity. 
The !99?'Tb(n, y) rate of this work is fitted as a function of temperature To (in unit of 10? K) with the standard form 
of REACLIB: 


Na < ov >(n,9)= exp(16.7928 + 0.0115T, ! — 1.18367; !/? +.4.0069T; ’? — 0.512773 + 0.031455? — 1.35771n Ty) (9) 


with the fitting errors are less than 1% in the range from Tg = 0.02 to Tg = 10. 


IV. ASTROPHYSICAL IMPLICATIONS 


In TP-AGB stars, the !?C pocket is formed following the third dredge-up (TDU) via !?C(p,y) ?N(8*v)1?C and 
generate neutrons by reaction !?C(o,n) that trigger the s-process [59-61]. The impact of !1??'Tb stellar neutron capture 
cross sections on main s-process nucleosynthesis was investigated using the stellar evolution code MESA [62]. The 
stellar structure evolution and the s-process nucleosynthesis are computed with the MESA star and net modules 
separately, hence they requiring less computing time and resources. A low mass star with initial mass 2M and 
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metallicity Z = 0.01 evolved from the pre-main sequence to near the end of the TP-AGB phase with nuclear network 
agb.net. In the !?C pocket formed after the 3rd TDU episodes, the temperature varies within 5.5-8.9x 10" K and the 
density increase from 96.0 g/cm? to 7.7x10? g/cm? during the whole interpulse phase. 

To treat the s-process, detailed nucleosynthesis calculations were performed with MESA net module. The network 
include 775 isotopes and 7344 reactions, and the temperature, density and light nuclei initial abundance are extracted 
from the !?C pocket of the precalculated stellar structure. The reaction rates corresponding to all proton, a, 8- 
decay, and neutron capture are provided by JINA REACLIB data tables, based on the 2017 updated version of the 
compilation by Cyburt et al. [63]. In addition, the initial element abundances are taken from Lodders et al. [64] for 
both structure and nucleosynthesis simulations. 
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FIG. 7: (color online) (Top) Time evolution of nuclear abundances on the test trajectory with T = 8.9x107 K and p = 7.7x10? g/cm’. 
(Middle) The abundance differences between the default case of the JINA REACLIB rate for 1599 Tb(n, 4)199 Tb, which is derived from 
KADoNiS-v0.3, and the case of ENDF/B-VIILO. (Bottom) Same physical quantities as in Middle one but the reaction rates of this work 
instead of ENDF/B-VIIL.0 data are used in calculations. 


The sensitivity of nuclear yields is evaluated to the reaction rate of °°Tb(n, 5)199Tb. Fig. 7(a) shows the time 
evolution of the abundances of nuclides around A ~ 160 along the s-process reaction path at constant temperature 
(T = 8.9x10" K) and density (p = 7.7x10? g/cm?) with the default case of the JINA REACLIB rate for °Tb(n, 
y)19? Tb, which is derived from the KADoNiS-v0.3 database. The abundance of !6?Dy, which has a half-live of about 
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50 days in stellar environment [65], decreases rapidly due to the considering 6-decay process to !9?Ho. The abundances 
of seed nuclei with initial values around A z 160 decreased significantly by neutron capture reactions before evolving 
to 7x10!? s because of the s-process flow from the dominant seed nuclei ?9Fe not yet overcomes the first s-process peak 
nuclei Sr and Ba, which act as the bottle necks in the reaction flow. On the other hand, the abundances of !®°Dy and 
164Fr increased obviously with the depletion of !99'Tb and !6?Ho from about 107—4x10? s. Then, the productions of 
all nuclides around A = 160 increase continuously up to 2-3 orders larger than initial values at the end of interpulse. 
Fig.7(b) shows the abundance changes between the KADoNiS-v0.3 and the ENDF/B-VIII.0 rate of °°Tb(n, y)!©°Tb 
used in the network calculations. In this case, the data of KADoNiS-v0.3 is smaller about by 40% than the one of 
ENDF/B-VIIILO. The abundances of 199:169'Tp, 160.161 Dy noticeably are affected by the uncertainty of cross section of 
the °°Tb(n, y)'©°Tb reaction. While, this propagation effects are less important for those nuclei heavier from !°°Tb. 
This trend is caused by achieving s-process reaction flow equilibriums [1, 66, 67] due to the high neutron fluence in 
the 8C pocket. Fig.7(c) shows the same physical quantities as in Fig.7(b) but with the !9?'Tb(n, 4)19"Tb rate of 
this work instead of the ENDF/B-VIII.0 data in calculations. As shown in the figure, the changes of abundances are 
much smaller than those of Fig.7(b) because of good agreement between two rates within 5%. 


V. SUMMARY 


We measured !?9?Tb(n, y) cross sections up to 1 MeV neutron energy by TOF method using the CgD¢ detectors 
at the CSNS Back-n facility. Resonance capture kernels are determined from 1 eV to 1.2 keV by analyzing the 
measured data with the R-matrix code SAMMY. The cross sections are obtained relative to the gold standard over 
the energy range of 2 keV to 1 MeV and agree well with the TALYS-1.9 theoretical calculations and the JEFF-3.3 
evaluation library within the experimental errors. On the other hand, the ENDF/B-VIILO evaluated cross section 
significantly overestimated to those of the present study. Maxwellian-averaged neutron capture cross sections (MACS) 
extracted from the present (n, y) data are in good agreement with the data of JEFF-3.3 and KADOoNiS-v0.3 [17] 
within 596 and 396-1196 of discrepancies in the whole energy range, respectively. However, the data of ENDF/B-VIII.0 
and KADoNiS-v1.0 [19] are significantly overestimated to the present experimental results by 40% and 2096. The 
sensitivity analysis of 9? Tb(n, y)'©°Tb reaction rate is investigated for the stellar evolution and nucleosynthesis of 
the 2M star model using the MESA code. Significant changes in abundances of 199:160'Dp, 160.161 ]yy are observed 
by comparing the stellar reaction rates of !9?Tb(n, y)!©°Tb between KADoNiS-v0.3, this work and ENDF/B-VIII.0 
evaluation. The present results show very small propagation to the more heavier elements for changing the rate of 
9Tb(n, y)'©°Tb because of the establishment of a reaction flow equilibrium between the s-process main component 
nuclei. 
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Appendix A. Resonance Parameters 


The resonance parameters obtained from fitting the measured capture yield of !??Tb(n, y) with R-matrix code 
SAMMY are compared with those from JEFF-3.3 evaluations in Tables VII and VIII. In the energy region from 100 
eV to 2 keV, the resonance parameters are not sufficiently reliable due to the worsening of the experimental resolution 
at Back-n facility. Therefore, only the resonance energies Eg and capture kernels k ( k = gV,Ty/(Tn 4- T4), g is 
the statistical factor) in this energy range are listed in VIII. Average resonance parameters, namely the average level 
spacing D;, the neutron strength function S; and the average radiation width < I, >), obtained in this work and the 
recommended data are listed in Table IX. The uncertainties in Tables of this work are from the fitting procedures. 
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Table VII Resonance parameters up to 100 eV neutron energy for ?'Tb(n, y) reaction. 


I^[meV] 


En [eV] 


Is[meV] 


En [eV] 
3.360 
4.990 
11.130 
14.500 
21.190 
24.520 
27.510 
33.810 
40.740 
43.630 
45.990 
50.070 
51.550 
53.990 
57.350 
58.630 
65.091 
66.484 
73.665 
76.318 
TE TOS 
78.598 
88.240 
90.310 


97.004 


JEFF-3.3 
L4[meV] Tg[meV] 
0.34 103.00 
0.08 103.00 
7.69 99.00 
0.19 105.00 
1.14 102.00 
5.32 116.00 
0.83 102.00 
2.61 98.00 
0.84 101.00 
5.90 97.00 
13.94 109.00 
1.91 96.00 
0.84 96.00 
0.83 78.00 
2.20 99.00 
1.59 96.00 
12.61 96.00 
3.52 98.00 
19.09 98.00 
6.92 108.00 
7.25 96.00 
2.68 85.00 
3.37 70.00 
6.84 90.00 
38.21 101.00 


0.71 


0.52 


1.59 


0.31 


3.48 


7.72 


1.17 


0.52 


0.31 


0.81 


0.98 


6.97 


1.27 


5.99 


2.44 


2.53 


0.97 


2.01 


3.97 


10.40 
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'Table VIII The resonance energies and kernels in this work are compared with the JEFF-3.3 
from 100 eV to 1.2 keV. 
'This work JEFF-3.3 
En [eV] k En [eV] k ER [eV] k En [eV] k 
109.24+0.22 8.24+0.34 111.56+0.27 3.04+0.217 | 108.98 8.19 111.26 2.76 
113.89+0.26 6.79+0.472 115.73+0.32 2.25+0.173 | 113.61 10.66 115.41 3.62 


118.83+0.28 2.45+0.225 128.46+0.25  0.63:-0.0465 | 119.06 | 5.68 128.13 0.53 


137.28+0.62 = 1.2840.267 138.25+0.45 2.960.366 137.43 1.45 137.89 3.57 


141.47+0.26 14.80+2.91 143.69+0.26 2.550.942 141.16 13.63 143.3 3.98 


152.90+0.24 9.670.747 155.95+0.28 6.080.479 152.52 7.88 155.43 13.48 


167.96+0.40 0.830.0827 170.13+0.30 3.500.258 | 167.47 0.77 169.63 3.1 


173.03+0.38 2.070.144 177.21+0.40 6.400.447 | 172.56 1.5 177.3 8.37 


186.00+0.53 0.470.043 197.78+0.44 6.730.414 185.42 0.8 197.36 7.45 


200.24+0.48 1.070.135 201.84+0.02 0.830.439 199.60 1.81 201.37 1.39 


211.15+0.54 0.73+0.0579 219.61+0.39 12.13+0.893 | 210.80 0.71 218.99 28.42 


226.33+0.78 0.970.109 229.15+1.60 0.410.102 225.90 0.79 228.4 0.38 


235.67+0.57 15.150.928 239.34+0.43 14.12+1.43 | 235.10 13.25 238.6 11.02 


242.17+0.45 5.410.516 245.28+0.16 0.15+0.0367 | 241.43 5.98 244.4 0.47 


251.93+0.61 1.90+0.18 254.75+0.58 0.79+0.0953 | 251.18 3.77 254.1 1.23 


263.24+0.71 6.890.519 269.05+0.61 5.16+0.329 | 262.90 10.68 268.39 4.15 


273.83+0.58 7.150.557 280.05+0.96 1.82+0.341 273.37 13.24 279.8 1.22 


281.86+0.69 6.61+0.803 284.77+0.53 17.03+3.04 281.3 6.71 284.14 20.96 


290.57+0.59 3.680.362 301.43+0.61 16.92+1.03 | 290.12 6.51 300.76 14.29 


306.22+0.65 10.400.813 312.88+0.65 14.52+0.968 | 305.49 9.36 312.11 17.66 


316.27+0.99 1.730.322 324.12+1.12 8.25+1.21 315.29 2.12 323.1 6.83 


326.58+1.26 10.4942.57 329.02+0.71 22.49+5.41 325.54 5.41 328.01 19.12 


332.86+0.64 21.9042.48 337.82+0.56 0.13+0.0832 | 331.94 241 336.92 1.08 


340.49+0.59 0.180.0908 346.16+2.18  12.3942.49 | 339.43 0.95 345.28 4.3 


348.50+1.58  10.6743.24 350.83+1.34 4.58+1.19 347.46 13.32 349.88 5.62 


continued on next page 
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Table VIII(continued) 


En [eV] 


358.904 


370.074 


375.374 


385.424 


409.994 


433.164 


442.324 


451.344 


455.154 


463.834 


475.374 


489.833 


499.214 


511.474 


521.884 


533.604 


547.174 


558.284 


566.044 


576.444 


581.014 


594.834 


600.414 


605.903 


615.814 


627.474 


639.304 


649.114 


-0.78 


-1.79 


E1:21 


-0.82 


-1.22 


eal 


2.07 


-4.06 


-2.46 


-0.02 


-1.47 


-0.87 


-3.03 


-1.65 


- 1.58 


- 0.88 


-2.92 


-4.99 


-2.01 


-2.19 


-2.24 


-2.48 


-2.16 


-3.19 


-0.94 


-3.03 


-3.02 


-1.49 


660.533 


-1.76 


This work 


k 


39.253-2.52 


6.75:-2.2 


13.5042.54 


2.830.345 


4.82+0.868 


6.32+0.691 


7.75+3.04 


5.19:-2.31 


8.02+3.1 


0.26£0.25 


22.85+2.47 


31.60+5.5 


8.99+2.87 


4.770.695 


25.19+2.91 


14.933-2.81 


5.3841.75 


4.802.44 


33.9145.72 


2.97+1.02 


5.56£1.32 


24.69+8.43 


17.77+45.24 


4.63£1.68 


27.312:9.28 


9.57+3.67 


11.98+2.34 


23.025.35 


4.010.684 


Er [eV] 


367.324 


372.774 


319.404 


403.204 


429.984 


440.264 


444.974 


453.394 


458.694 


466.154 


482.734 


495.584 


504.154 


518.434 


529.194 


544.844 


553.364 


561.604 


571.473 


518.4243 


592.314 


597.474 


603.104 


610.504 


619.304 


630.284 


644.794 


658.224 


-1.40 


-0.88 


-0.82 


-1.61 


-2.88 


-2.65 


-2.60 


-1.91 


-1.06 


-2.81 


-2.17 


-2.18 


-1.35 


-3.54 


-1.47 


-2.70 


-2.92 


-3.53 


-2.78 


-3.35 


-3.82 


-3.25 


-2.68 


E3.57 


-2.74 


-3.80 


-3.11 


-3.85 


663.244 


-1.33 


k 


10.58+1.54 


22.02+3.26 


18.89+1.62 


8.001-1.55 


0.570.218 


10.52+2.66 


5.814-1.58 


10.79+3.73 


9.02+1.2 


2.19-£0.846 


3.49+4.38 


10.62+2.69 


13.48+1.6 


5.2341.87 


17.40+2.07 


4.07+1.42 


8.911.85 


12.98+4.92 


10.743-2.55 


4.06:-1.49 


23.21+6.52 


16.48+6.58 


9.54+3.21 


16.60+3.29 


7.3341.91 


17.25+4.37 


11.45+3.57 


1.070.358 


4.110.636 


Er [eV] 


358.70 


368.64 


374.14 


384.80 


409.02 


432.24 


442.60 


450.98 


453.10 


463.32 


474.00 


488.52 


497.77 


510.06 


520.78 


532.83 


546.32 


557.02 


564.71 


575.70 


579.68 


593.60 


600.02 


605.50 


615.34 


627.03 


637.60 


648.35 


659.92 


JEFF-3.3 
k Er [eV] 
29.9 366.32 
3.67 371.81 
11.13 378.16 
3.4 404.5 
6.25 429.1 
6.29 439.96 
9.39 444.45 
1.92 452.8 
3.7 457.7 
113 4067.87 
12.53 482.07 
37.82 493.97 
6.55 502.96 
12.95 517.37 
21.44 527.79 
18.45 544.5 
1.9 551.8 
Cl 560.47 
16.34 569.73 
4.59 577.39 
14.88 592.39 
35.43 597.52 
11.05 602.71 
2.31 608.12 
50.31 620.49 
22.74 630.29 
8.15 644.28 
31.82 654.5 
8.03 662.6 


k 


6.69 


16.48 


20.49 


18.3 


0.64 


11.49 


3.42 


15.47 


13.7 


33.72 


27.77 


12.96 


14.56 


3.68 


19.16 


9.64 


8.88 


8.11 


12.92 


2.77 


17.81 


19.04 


13.72 


8.64 


11.09 


8.11 


7.86 


2.46 


5.29 


continued on next page 
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Table VIII(continued) 
This work JEFF-3.3 
En [eV] k En [eV] k En [eV] k | Em [eV] k 


678.45+2.51 4.75+1.06 681.88+2.86 7.4142.1 677.70 4.55 682.51 16.06 


684.57+2.23 12.8642.91 687.71+2.97 18.41+4.44 | 683.60 34.2 686.64 14.48 


694.55+4.29 8.902-2.98 702.00+3.12 15.57+3.21 693.48 7.61 700.57 14.2 


708.96+2.52 12:70£2.27 720.87+1.91 8.452.84 707.47 6.79 719.68 33.22 


723.53+0.40 3.99+2.51 127.88+4.88 3.51£1.75 723.08 1.94 726.97 6.11 


732.89+2.91 27.50+4.62 740.27+5.74 6.82+2.59 731.75 25.83 738.5 9.01 


748.01+3.48 4.48+1.44 750.91+3.58 4.70+1.61 747.32 7.29 751.46 10.19 


794.69+3.68 5.9141.77 767.45+2.14 57.03£7.81 | 753.90 5.15 765.2 40.77 


768.77+6.80 25.421-22.3 783.42+4.63 6.0842.81 769.80 19.03 782.5 3.92 


787.44+2.68 37.048.35 790.76+2.32 21.69+5.76 | 787.70 28.37 793 6.31 


802.50+4.95 12.413.37 810.17+5.36 10.99+5.99 | 802.00 23.45 810.2 16.41 


814.44+7.91 15.656.87 824.44+2.35 35.91+3.31 815.70 5.22 823.9 22.62 


837.62+4.69 2.010.792 845.013.69 39.72+7.69 | 837.00 3.88 845.3 25.67 


849.34+2.22 31.803-6.76 853.81--3.24 19.32+5.49 | 850.50 26.82 855.9 6.53 


859.26+4.33 8.40+3.06 872.3319.42 3.2442.78 861.40 8.88 869.8 4.12 


876.21+6.60 14.09+5.05 884.20+4.27 25.68+7.24 | 875.40 5.39 884.7 21.25 


896.76+9.85 6.214:5.75 901.552-3.01 25.354 


7 896.40 9.01 901.5 24.29 


911.67+10.04 7.09+4.4 925.28+2.25 30.40+5.38 | 921.50 4.86 924.9 32.15 


935.80+2.14 138.65+11.8 951.33+3.99 5.43+1.82 936.00 42.09 951.3 17.54 


955.25+6.16 2.95+1.56 966.97+2.56 21.96+6.41 | 955.30 6.37 968.4 10.9 


971.69+2.67 21.47+4.59 979.62+3.80 29.03+6.75 | 972.60 17.65 980.2 29.62 


994.59--2.90 42.17+6.68 998.76+5.29 24.14+1.07 | 995.20 43.79 1000.7 20.79 


1005.19+78.01 1.3640.14 1015.76+77.30 6.810.68 1005.00 5.74 1016.2 12.62 


1028.42+78.92 1.0341.03 1036.81+78.88 0.280.283 | 1028.00 6.57 1036.3 17.27 


1040.663-79.00 | 0.043-0.0372 1050.39+79.47 0.06+0.0572 | 1040.20 23.26 1050 33.11 


1057.37+79.81 0.030.031 1068.11+79.95 4.79+4.76 | 1057.00 31.37 1067.7 14.64 


1097.87+81.84 0.170.173 1105.74+82.28 0.310.313 | 1097.60 24.61 1105.5 17.37 


1124.31+83.28 0.930.929  1132.416-83.67  0.77£0.77 1124.10 20.15 1132.2 21.64 


1142.814-85.02 2.980.298 1148.20+84.49 0.890.888 | 1142.70 7.42 1148 21.73 


continued on next page 
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Table VIII(continued) 
This work JEFF-3.3 
Er [eV] k En [eV] k En [eV] k En [eV] k 


1157.19+85.48 10.46410.5 1172.964135.74 39.86+46.8 


1157.00 13.28 1172.8 51.88 


1184.70+87.74 4.2324.21 1192.36486.92 1.14+1.14 


1184.40 17.68 1192.2 15.7 


Table IX Average resonance parameters of this work and the recommended data from 


S.F.Mughabghab [51] and ENDF/B-VIII.0 [52]. 


Parameters ThisWork ENDF/B-VIILO — S.F.Mughabghab 
So (107^) 1.53+0.02 1.207 1.55+0.15 

Sı (10~*) 1.83+0.05 1.480 1.86+0.09 
S2(10~*) 1.55+0.21 1.090 1.41+0.24 
<T, >o(eV)  0.09840.002 0.097 0.101+0.002 
<T, >ı(eV) 0.055+0.003 0.097 0.050+0.004 
<T, >2(eV) 0.098 0.097 - 
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